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Interaction of domain walls (DWs) in ferromagnetic stripes is studied with relevance to the forma- 
tion of stable complexes of many domains. Two DW system is described with the Landau-Lifshitz- 
Gilbert equation including regimes of narrow and wide stripes which correspond the presence of 
transverse and vortex DWs. The DWs of both kinds are characterized with their chiralities (the 
direction of the magnetization rotation in the stripe plane) and polarities (the magnetization ori- 
entation in the center of a vortex and/or halfvortices), hence, their interactions are analyzed with 
dependence on these properties. In particular, pairs of the DWs of opposite or like chiralities and 
polarities are investigated as well as pairs of opposite (like) chiralities and of like (opposite) polar- 
ities. Conditions of the creation of stationary magnetic bubbles built of two interacting DWs are 
formulated with relevance to the situations of presence and absence of the external magnetic field. 
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I. INTRODUCTION 

Recent growth of the interest in complexes of domain 
walls (DWs) in quasi- ID ferromagnets is due to hopes for 
miniaturization of information registers and logical de- 
vices based on magnetic nanowires, [lHj], and nanorings 
@. Currently, the main technological effort is focused 
on the (straight or curved) magnetic nanostripes whose 
spin structure is more complex than described within the 
ID idealization. It is because nanostructures of the best 
quality are obtained in this form using lithography meth- 
ods [6|. Mutual interactions of the DWs influence the 
stability of a record of bits encoded in a string of the 
magnetic domains Q. 

In the present paper, I study interactions of DWs in 
ferromagnetic stripes of nonzero thickness including two 
regimes of the stripe widths; narrow stripes correspond 
to the so called transverse DWs while wider stripes to 
the so called vortex DWs 041] ■ It is shown that both 
kinds of the DWs are exact stationary solutions to the 
Landau-Lifshitz-Gilbert (LLG) equation in 2D (different 
realizations of the cross-tie DW 0[H|]), however, their 
orientations in the stripe plane are distinguished by dif- 
ferent boundary conditions. Both transverse and vortex 
DWs are characterized by two features; a chirality (the 
direction of magnetization rotation in the stripe plane) 
which is clockwise or counterclockwise one and a polar- 
ity (the direction of magnetization alignment in a vor- 
tex/halfvortex center) which is positive or negative one 
[H G3]. Pairs of the DWs of the opposite (like) both 
the chiralities and polarities are found to constitute ex- 
act stationary solutions to LLG in the absence of external 
field, hence, they do not interact. The interactions of the 
DWs of opposite chiralities and of like polarities (as well 
as of like chiralities and of opposite polarities) are stud- 
ied within a perturbation calculus previously developed 



with relevance to DWs in ID ferromagnetic wire and in 
critical systems |l4(. Evaluating the dependence of en- 
ergy of the two-DW systems on the distance between the 
DWs, I analyze the creation of stationary bound states 
of two DWs (magnetic bubbles). 

In the presence of an external magnetic field, even DWs 
of the opposite chiralities and polarities (of like chirali- 
ties and polarities) interact. The field- induced collision 
of such DWs is studied applying a method developed to 
ID systems in [15]. For this aim, an extension of the 
(dissipative) evolution equation of magnetization is per- 
formed in a way to describe the dynamics in the limits of 
large positive and large negative values of time. The col- 
lision is found to be accompanied by the reflection of the 
DWs, hence, by the creation of stable magnetic bubbles. 

In section II, the single DW solutions to the LLG equa- 
tion in magnetic stripes are analyzed. The interactions of 
the transverse and vortex DWs of opposite chiralities and 
of like polarities (as well as of like chiralities and of oppo- 
site polarities) are investigated within the perturbation 
calculus in section III while the field-induced collisions of 
the DWs of opposite (like) chiralities and polarities are 
studied in section IV. Main conclusions are summarized 
in section V. 



II. SINGLE DW IN FERROMAGNETIC STRIPE 

I consider stationary DW solutions to the LLG equa- 
tion in 2D 
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The first term on the r.h.s. of ([T]) relates to the ex- 
change interactions while the second (Zeeman) term de- 
pends on the external magnetic field H = (H x ,Q,0), 
thus, 7 denotes the gyromagnetic factor. The constant 
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/3i determines strength of the easy-axis anisotropy and 
i = (1,0,0) indicates the long axis of the stripe. In 
nanomagnets with structural disorder (permalloys are 
the most popular materials for nanostripcs), the effec- 
tive spin anisotropy is dominated by dipolar interactions, 
hence, it is a shape dependent and mainly surface ef- 
fect. Then, for stripes of big enough cross section, the 
bulk anisotropy is negligible, j3\ <C J/w 2 , where w de- 
notes the stripe width. Since (fTJ is valid only when the 
constraint |m| = M is satisfied, one writes equations of 
the unconstrained dynamics equivalent to ([T]). Introduc- 
ing m± = m y ± im z , one represents the magnetization 
components using a pair of complex functions g(x,z,t), 
f(x,z,t). The relation between the primary and sec- 
ondary dynamical variables 
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ensures that |m| = M. Inserting @ into (flj leads, fol- 
lowing the Hirota method for solving nonlinear equations 

[Bill, to 

[-iD t + J(D 2 x + Dl)+aD t ]r-g 

-( 1 H x + p 1 )f*g = 0, 
[-iD t - J{Dl + Dl)+aD t ]t-g 

+(- 7 ff x + /?!)/*<? = 0,(3) 
(Dl+D 2 z )g.g = 0, (D 2 X + D 2 z )f* ■ f* = 0, 

where D t , D x , D z denote Hirota operators of differenti- 
ation 

D™b(x, z, t) ■ c{x, z, t) = 

(d/ dx - d/ dx') n b(x, z, t)c{x ,z',t )\ x=x , >z=z , t=t , . 

For H = 0, stationary single-DW solutions to (J3|) take 
the form 
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FIG. 1. DW configurations; a) a transverse DW, b) a vor- 
tex DW. In the upper draws, arrows indicate magnetization 
alignment. 
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and Kek ^= 0. I denote k = k + ik , q = q + iq . 
Assuming one of the DW ends to be centered at x = 0, 
(then u = e 1¥> ), the relevant magnetization profile [the 
single-DW solution to (fTJ] is written explicitly with 

m+{x, z) = Me i(ip+k " x+q " z) sech[k'x + q z), 
m x (x, z) — — Mtanh[fc x + q z]. (6) 

I assume q ^ since, in the case q =0, the DW 
states are similar to DWs in ID and cannot exist in 
absence of the bulk anisotropy [H], HH, G3- Defin- 
ing 9 = arctan(g jk ), following <j5j) , one finds k — 



-q tan(0) and k 
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/3i/{J[l + tan 2 (6»)]}. Also, I 



assume the magnetization alignments on both the stripe 



edges to be similar, the ordering is symmetric with re- 
spect to the line z = w/2. It leads to the condition 
k (—wq /k ) + q w = mr where n = 1, 2, . . . and, finally, 
to q" = mr/{w[l + tan 2 (6»)]}. 

Additional boundary condition is related to minimiza- 
tion of the surface (magnetostatic) energy and it discrim- 
inates between different values of (p, n, and 8. I evaluate 
the dependence of the energy of the DW on these param- 
eters using the Hamiltonian % — T-Lq + TLz, where 
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T-Lz = -7 H • m, 
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(Hz denotes its Zeeman part). Total energy of the DW 
E = Eq + Ez + Eb is the sum of the bulk energy 
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E a + E z , defined by E ( Z) = j!!^ J™ H (z)dzdx, and 
of the boundary energy Eg. The contribution Eb is de- 
termined referring to a theorem by Carbou who proved 
that the magnetostatic energy of any ferromagnetic el- 
ement of finite thickness r; 1/X 2 J s (m ■ n) 2 ds tends to 

I/A2 J as ( m • n') 2 d/ with t — y 0, [ll]. Here, S denotes 
the surface of the bulk ferromagnet and dS denotes the 
boundary of the base of its solid, n is normal to the mag- 
net surface, n denotes the unitary vector outward to the 
line of the base boundary. The coefficient A2 scales with r 
and with a diameter w following A 2 ~ A 2 /[r| log(r/W)|], 
(with relevance to a stripe, w represents its width j8l.ll9|). 
Writing the formula for Eb, one has to notice that, the 
Carbou theorem is not strictly applicable to systems with 
open boundaries, (infinite stripes), however, it shows 
some tendency in ordering at the stripe edges. In par- 
ticular, it indicates that the magnetostatic interactions 
in flat magnets induce more than one hard directions of 
magnetization parallel to the main plane. I propose to 
effectively describe the magnetostatic energy in the form 
of an integral over the stripe edge 
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where, in contrast to [19j, the effect of anisotropy relat- 
ing to the long axis of the stripe is included via the term 
depending on Ai. By analogy to systems with bound- 
aries of finite length that satisfy the Carbou theorem, 
this coefficient is expected to scale with the stripe width 
following Ai ~ A 2 jw. More detailed estimation of Eb is 
performed in Appendix A, (see also (2fjjV 

Inserting ^ into the Hamiltonian 0, one arrives at 



E (6,n) = 2JMw\j y[l + tan 2 (6>) 
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Estimating Eb, I divide it into two parts Eb = Ebi 
E B 2-, 
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The minimization of EB(<p,0,n) leads to ip — 0,7r inde- 
pendently of other parameters of the DW. It corresponds 
to the presence of halfvortices (halfantivortices) at the 
stripe edges (as shown in Fig. 1) . In the regime of narrow 
stripes w/t ~ 1, A2 <C Ai, (I/A2 is small, thus, Eq domi- 
nates over Eb), via minimization of Eq{6, n) the smallest 
possible value of n is preferable while the minimization of 
both the energy components Eo(0,n), Eb(0,8,ti) leads 
to the preferred values of the angle 8 — 0,n. Increase of 
the stripe width w with fixed thickness r results in the 
transition to the regime A2 > Ai. Furthermore, as men- 
tioned previously, for soft-magnetic alloys, Pi is negligi- 
bly small, Pi < J/w 2 , thus, from ©, E (8,n) = E (n) 
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FIG. 2. Scheme of the magnetization layout in double-DW 
configurations; a) pair od transverse (parallel vortex) DWs 
of opposite chiralities and of like polarities, b) pair od trans- 
verse (parallel vortex) DWs of like chiralities and of opposite 
polarities. 



is independent of w, and Eb becomes comparable to 
E Q . For this case, via the minimization of Ebi{0,7i), 
the biggest possible value of tan 2 (6>) and the smallest 
value of n are preferable, whereas, for 8 ^ 0,tt, the min- 
imization of EB2(0,9,n) indicates big values of n to be 
preferable. Hence, the transition from the DW state of 
n = 1, 9 = 0,7r to a state of n = 2, 8 ^ takes place 
with increase of w. For n = 2, the condition tan 2 (6>) > 1 
has to be satisfied in order to £^(0, 8, 2) < ^(0, 0, 1). 

The state n = 1 and 8 — 0, tt corresp onds to k = q = 
0, \q I = tt/w, \k I = \Jtx 2 jw 1 + Pi/ J and it is called a 
transverse DW. With relevance to the state n — 2,1 take 
Pi = 0, (the shape anisotropy does not affect the dy- 
namical equation while it is completely included in the 
boundary condition, the minimization of Eb), and \8\ to 
be close to its infimum \8\ w 7r/4. The resulting mag- 
netization structure corresponds to |q | ~ |<7 |«|fc|w 
\k « tt/w and one calls it a vortex DW. The polarity of 
vortex (transverse) DW, (the magnetization orientation 
in the center of vortex/halfvortex, parallel or antiparallel 
to y axis), is determined by value of tp while sign(g") de- 
termines its chirality, (the magnetization rotation in the 
stripe plane, clockwise or anticlockwise to y axis). 



III. INTERACTION OF STATIC DOMAIN 
WALLS 

Let me analyze the interaction of transverse (vortex) 
DWs within a perturbation calculus developed earlier 
with relevance to the interactions of DWs in ID ferro- 
magnet [3]. In the present section, I focus on the sys- 
tems whose interacting DWs are of opposite chiralities 
and of like polarities (or of like chiralities and of opposite 
polarities) visualized in Fig. 2. I refer to the remaining 
cases (the DWs of opposite chiralities and polarities as 
well as the DWs of like chiralities and polarities) in the 
next section. 

Locally, in the vicinity of the center of jth. DW, (j = 
1,2), one can write the magnetization in the form 

m(x, z,0) = m^(a;,2) + Sm^(x, z), (11) 
where mW denotes the stationary single-DW solution to 
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(x, z) = Me i ^- Vj+5j ^ k ( x - x °^+i 2 1} 

xsech{(ij[A: (x — xqj) + q z]}, 



mS^^x, z) = — MtanhjcTj [k (x — XQj) + q z]} (12) 

while Sm^ denotes a perturbation due the presence of 
another DW. Here Si = S 2 = ±1, \aj\ — 1, (J\ = — cr 2 , 
and, for k > 0, a } ■ — 1 corresponds to DW of the head- 
to-head type while crj = — 1 to the tail-to-tail type. The 
phases satisfy tpi — ip2 or <pi = ip2 ± tt. I refer to the 
case Si = 82 and y>i = <f2 as to the pair of walls of op- 
posite chiralities and of like polarities while to the case 
81 = S 2 and ipi = if2 i 7r as to the walls of like chirali- 
ties and of opposite polarities. In terms of vortex DWs, 
I focus my attention to the interaction of mutually par- 
allel vortex DWs while taking the single DWs profile in 
the form ([12")) excludes its use to non-parallel DWs. The 
non-parallel vortex walls overlap at one edge of the stripe 
much stronger than at the second edge, thus, an effective 
description of their interaction can be performed with ID 
model of Ref. [l4[ applied to the stripe-edge magnetiza- 
tion. In order to fulfill the constraint |m| = M, I take 
the perturbation of the form 
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The plus and minus relate to the bubble magnetization 
parallel and antiparallel to x-axis, respectively. I stress 
that ([14")) is irrelevant to the case of DWs of opposite (like) 
chiralities and polarities, Si = —S 2 , since, for this case, 
I m I = M was not satisfied. The perturbation of states 
of the interacting DWs should be small. According to 
(TIB")) , this condition demands the interacting DWs to be 
far-enough from each other, (functions of the r.h.s. of 

(US)) satisfy | ± m ( x k) /M - 1|, \m ( y } /M\, \m (k) /M\ < 1 if 
the distance between walls is larger than their width). 

With relevance to the transverse DWs, insertion of 
(fT4")) . with \q'\ = \k"\ = 0, into the Hamiltonian © 
and integration over the stripe area {x G (—00,00), 
z £ [0,ui]} leads to the following dependence of the en- 
ergy of a DW pair on the distance of separation of the 
walls 
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where a = k (xq2 — xqi), and / ± (a, 9) represent integrals 
written in Appendix B. The upper signs correspond to 
the pair of transverse DWs of opposite chiralities and of 
like polarities, the case ip± — (f2, while the lower signs 
correspond to the pair of transverse DWs of like chirali- 
ties and of opposite polarities, the case ipi = (p2 ± tt. Up 
to a constant of proportionality, the above energy depen- 
dences on a are similar to the ones of the pairs of DWs (of 
opposite chiralities and of like chiralities, respectively) in 
ID ferromagnet [14j . 

For parallel vortex DWs, inserting ([i"4)) . with \q ' )| = 
\k ' -*| = n/w and /3i = 0, into the Hamiltonian and 
integrating it over the stripe area, one arrives at 

E (a) = MwJk 7 ± (a, 1) 

= 7rMJ-csch J (a){[36acosh(a) — 12acosh(3a) 
6 

— 24sinh(a) — 5sinh(3a) + 3sinh(5a)] 
=F24cos(a)sinh 2 (a)[-2a + sinh(2a)] (16) 
— cos(2a)[24acosh(a) — 18sinh(a) — 2sinh(3a)]}, 

where upper signs correspond to the pair of DWs of op- 
posite chiralities and of like polarities while lower signs 
to the pair od DWs of like chiralities and of opposite 
polarities. The energy dependences on the DW sepa- 
ration distance (|T5l) . (IT6l) are plotted in Figs. 3a and 
3b, respectively. The character of extremum (minimum 
or maximum) of the function Eq(o) at a = indicates 
that the interaction of the (transverse or vortex) DWs of 
opposite chiralities and like polarities is attractive when 
their separation distance is short while the interaction of 
DWs of like chiralities and opposite polarities is repul- 
sive in this case. To be precise, here, I call the distance 
between DWs short when it is close to their width since 
the perturbation calculus is applicable to \a\ > 1. 

In the case of presence of an external longitudinal field 
H x ^ 0, the energy dependence on the distance of DW 
separation E (a) + Ez(a) deviates from these in Figs. 3a 
and 3b due to non-zero Zeeman term, however, such a de- 
viation is small for a close to zero, similar to one in Fig. 
1 of Ref. [3]. Moreover, the Zeeman term is arbitrarily 
small when restrict considerations to weak-enough exter- 
nal field. With this restriction, due to attractive interac- 
tion, an (field-induced) collision of transverse or vortex 
DWs of opposite chiralities and of like polarities is ex- 
pected to result in their mutual annihilation, according 
to simulations for the transverse DWs [U ■ Corre- 
spondingly, due to repulsive interaction, the collision of 
transverse or vortex DWs of like chiralities and of oppo- 
site polarities leads to their mutual reflection. Both the 
DW annihilation and nucleation have been observed in 
ferromagnetic nanorigs [23|, [24[ ■ 

Unlike in the case of transverse DWs, the interaction 
of vortex DWs changes its character from attractive to 
repulsive (or vice versa) with increasing the distance be- 
tween the DWs. 
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FIG. 3. Energy of a pair of DWs with dependence on the 
distance of their separation; a) transverse DWs, b) vortex 
DWs. Solid line - DWs of opposite chiralities and of like 
polarities, dotted line - DWs of like chiralities and of opposite 
polarities. 



IV. FIELD-INDUCED COLLISION OF DOMAIN 
WALLS 

In the present section, the interaction of DWs of oppo- 
site chiralities and polarities (as well as of like chiralities 
and polarities) is analyzed. It is shown that such pairs 
of DWs create static bubbles of magnetization which are 
exact stationary solutions to the LLG equation. Hence, 
the walls do not interact in absence of any external field. 
The application of a magnetic field in the opposite direc- 
tion to the bubble magnetization induces motion of both 
the DWs towards each other and, eventually, their colli- 
sion. It is accompanied by an interaction of the DWs due 
to a dynamical deformation of them which cannot be de- 
scribed in the framework of the above used perturbation 
calculus. 

Stationary two-DW solutions to ([T]) have been found 
in the absence of anisotropy, for /3i — 0. Under this 
condition, for H = 0, inserting the ansatz / = 1, g = 

ue k 1 x+q 1 y ± ue k 2 x+q 2 y ^ Qne findg ^ re l ations 

fc l(2) + 9l(2) = °J fc l fc 2 + <7l<?2 = 0. (17) 

Let kj = kj' + ik/', qj = q/ + iq/'. For double- DW so- 
lutions, k x = — k 2 = ±ir/w, the above conditions lead 
to q[ = -ql, \q' m \ = n/w, k'[ {2) = q' m = for 

transverse DWs, and to k Y k 2 — k 1 k 2 , q±q 2 = q±q 2 , 

1^1(2)1 = 1^1(2)1 = 1^1(2)1 = 11 1 W ^ 0r vor t ex DWs. 

Concerning nonstationary DW solutions to the LLG 
equation for H / and a / 0, I notice that solitary- 
wave solutions to the LLG equation are relevant in the 
limit of large positive values of time only, whenever mag- 
netostatic effects which lead to the Walker breakdown ( an 
easy-plane anisotropy in the case of ID ferromagnet) are 
included [25T - [27| . Hence, studying the dynamics of a pair 
of DWs which are stationary at the initial moment, one 
cannot avoid an effect of non-adiabatic switching on of 
the magnetic field. On the other hand, standard ap- 
proach to the problem of soliton collision is based on the 
analysis of asymptotics of two soliton solutions in the 
limits t — > ±oo, e.g. it has been used with relevance to 



collisions of spontaneously propagating topological soli- 
tons (DWs) in ID ferromagnet (in absence of external 
field and dissipation) [l|| . In the presence of dissipation 
this method fails since solutions to dissipative equations 
of motion become unphysical in the limit of large nega- 
tive values of time, because of divergence of energy in this 
limit. In particular, Zeeman energy of any ferromagnetic 
stripe diverges with t — > — oo whereas domains magne- 
tized parallel (antiparallel) to the magnetic field grow 
(diminish) with time. In view of the purpose of studying 
field-induced collisions, these facts motivate extension of 
the dynamical system (J3J within a formalism applicable 
to the limits of large positive and negative values of time. 

The method of extension of the dynamical system has 
been developed in Ref. 15]. It is based on a Bateman's 
idea of doubling the number of freedom degrees when 
including the dissipation into the standard Lagrange for- 
malism (with relevance to damped harmonic oscillator) 
[28| . I change secondary dynamical equations © replac- 
ing g, g* , f, f* with novel fields of the corresponding set 
gi, g 2 , f 2 , /* and of the set of their c.c.s. in a way that 
the resulting extended system of the equations of motion 

aD t f" ■ gi - jH x f*gi = iD t f* ■ gi, 

{D 2 X + Dl)f{ ■ 9l = 0, (18) 
{D 2 X + D 2 z ) gi .g 1 = 0, {Dl + Dl)ft ■ n = 0- 

and of their conjugates [that differ from (fT8"|) by the sign 
of the dissipation constant a] 

- aD t f 2 ■ 92 - lH x f 2 g 2 = iD t f 2 ■ g 2 , 

(Dl + D 2 z )f;-g 2 = 0, (19) 
(Dl + Dl)g 2 ■ g 2 = 0, (D 2 X + D 2 z )f* ■ f * = 

is, in a formal sense, symmetric with respect to the rever- 
sal of the arrow of time, (although, dissipative dynamics 
of any physical system is irreversible). In the system 
of eight equations; (fT5|) - (fTi?)) and of their c.c.s, gi(fi) is 
not a c.c. to g 2 (f 2 ), while, comparing (|T5)) and (|19l) . 
one sees that g 2 (x, z, t) [f 2 (x, z, t)] can be obtained from 
gi(x, z, t) [/i(x, 2, t)) via changing the sign of its param- 
eter a. Upon the change t — > —t, the system of the 
novel equations transforms into itself if one accompanies 
this operation by the transform of the novel dynamical 
variables g^ 2 ) -> /2(i), /i(2) -> -52(1)- The above trick 
is an analogous construction to classical and quantum 
formalisms for description of dissipative systems in the 
whole length of the time axis despite the divergence of 
excitation energy in the limit of large negative values of 
time (non-equilibrium Green functions, thermo-field dy- 
namics, rigged Hilbert space) (29[, which are all based on 
the concept of Batcman. The equations (TTHJ) and their 
c.c.s., which determine the magnetization dynamics for 
large positive values of time (in particular, for t — > 00), 
contain the differentials of the functions g±, g\, f\, f*. 
Therefore, the magnetization vector should be expressed 
with these functions in the relevant time regime. Writing 
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the magnetization in the form 

m + = 7w 1 * 77 » m * = M 771 ; — rrr ( 20 

ensures that their components satisfy |m| = M, m x = 
to*, and they reproduce ([2]) for a = 0. In the limit t — >• 
— oo, one can analyze the evolution of the magnetization 
with vector m defined as 



2M 



TO + = 



fZ/92 + P2//2 



-M|/^4^(21) 
/ 2 /52 + S2//2 



The single- DW solution to ([T8 ]) -([T9 ]l is of the form 
f 2 = e~'*/ 2 , gx — ue kx +i z - lt / 2 1 which leads to the mag- 
netization profile 

m+(x,2,t) = Me i l« ,+i "( I - IO ' +9 " 2 - | " i ' 

xsech[fc (a: — xo) + q Z — 1 1], 
m x (x, z, t) — — Aftanh[fc (x — xq) + q z — I t], (22) 

and represents the translationally and rotationally mov- 
ing transverse DW of \k \ — \q | = ir/w, \k \ = \q \ = 
or vortex DW of \k ^ -* = |<? )| = ir/w. Here, Z = 
Rel = jH x a/(l + a 2 ), l" ee Iml = jH x /(l + a 2 ). 
The nonzero I corresponds to the DW evolution in 
strong enough magnetic field that exceeds significantly 
the Walker-breakdown value \H X \ 3> Hw, [3(1 IMj]. Be- 
low the Walker breakdown, in a weak external field, the 
rotation of the magnetization about x axis is suppressed 
by the minimization of surface energy, (a strong effec- 
tive bi-axial anisotropy 32]). The relevant dynamics is 
described with a differing from ([T]) primary evolution 
equation and differing from (fT5)) - (flTj| secondary equa- 
tions which are obtained by exchanging the l.h.s. of ([1]) 
and r.h.s.s. of p73|) - (ITOl) into zero [lfj. In this weak- field 
regime, I = jHx/a, I = 0. In the intermediate re- 
gion, above the Walker breakdown \H X \ > Hw, the field- 
induced motion of the DW is accompanied by a dynamic 
transformation of its structure [33| . The DW oscillatory 
changes between the transverse and vortex ones 34j . My 
further considerations focus on the regime \H X \ < Hw, 
then the motion of a single DW is a simple translation. 

For H ^ 0, two-DW solution to modified (by taking 
their l.h.s.s. equal to zero) Eqs. (|T5]) - (| 1 9 1) takes the form 



9l =u(e klX+qiZ ±e k 



k- 2 x+q 2 z^-~fH :c t/ (2a) 



= e -7#*t/(2a) 



(23) 



92 =u(c 



_j_ e fc 2 x+g 2 z^g7-Wxt/(2a) 



where plus (minus) corresponds to states of a pair of the 
DWs of like (opposite) chiralities and polarities, while the 
parameters fci(2), <?i( 2 ) satisfy the conditions (fl~7|) . Hence, 
for k 1 — —k 2 and q l = —q 2 , (parallel DWs); k l = —k 2 
and q 1 = —q 2 ■ It should be emphasized that the above 
single-DW and double-DW solutions to (jT5)l satisfy the 
original system ([3]). 



Let 

T)j(x, z,t) ee (jjTr/w(x — xqj + 9z) — jH x t/a, 
fjj(x, z, t) ee <7jir/w(x — xqj + 9z) + jH x t/a, (24) 
£j(x, z) ee 5jir/w[6(x — Xqj) — z], 

with ax = —o 2 = 1, 5i = —62, 1^1(2)1 = 1) and 9 = for 
transverse DWs while 9 = 1 for vortex DWs. For H x > 
and rjk <C to ~ 0, here k ^ j, (k,j = 1,2), we find the 
distant-future limit of the magnetization (|2U)) 



lim 

t— ► 00 



to+ = mf ee (±1) J - 1 2M 



ve Vj gifj 

1 + e 2 % 



lim to 2 

t— >oc 



m 



U) 



=2>7j 



M- 



1 



3 2% 



(25) 



Identifying the parameters xoj with the DW-center po- 
sitions x i = — X02, [In (it) oc <Ti + i<5i# = — <r 2 — iS 2 9], I 
introduce the restriction onu; v = 1 or « = — 1. The 
magnetization profiles (|25l) describe the motion of well 
separated DWs of the type presented with (|22"|) and corre- 
spond to the limit t — >• cxd of the DW solutions to the LLG 
equation in 2D for \H X \ < Hw by Slonczewski (27ll35j. 

In the distant-past limit, I describe the magnetization 
evolution using the field m. Following ff2T|) . for fjj -C 
fj k ~ and j ^ k, 



lim rhj 

t— >— 00 



- 0) 

mi' 



(±1 



,fe-i 



2M 



1 + c 2f i<* ' 



lim to x = m)P — —M 



1 - e 2r >* 
1 + e 2 ^ 



(26) 



Considering the collision of DWs which are infinitely dis- 
tant from each other at the beginning of their evolution, 
I determine the magnetization dynamic in the limit t — > 
—00. via inverting the propagation direction of the kinks 
of m and reversing the arrow's head of the field vector rh. 
Utilizing the properties m^P(x + xat, z,0)e 1 ^ k ^ x '~ z " 1 = 
m ( i\-x + a^-^Oy^*' 2 ), m { x j) {x + x okl z 7 0) = 

-mi\—x + x ok ,-z,0), for r)j > r/ k - 0, (j k), I 
arrive at 

m+(x,z,t) = -mf ] (-x + 2x Qk , -z,t)e i2 ^ x ' z '> ,(27a) 
m x (x, z, t) = —m:p (—X + 2xQ k , —z, t). (27b) 

The applicability of the above procedure to the asymp- 
totic evolution of a single DW can be verified noticing 
that any single-DW solution to (fl"g ]) -(fig j) . [Eq. ([52"])]. sat- 
isfies 



m+(x, z, t) = -fh+(—x + 2x , -z, t)e l2 ^ x - z ^ 
m x (x, z, t) = fh x (-x + 2x , -z, t). 



(28a) 
(28b) 



According to (j25l) . (|27a[) - (I27b|) . two initially closing up 
DWs of opposite (like) chiralities and polarities have to 
diverge after the collision. The one of the colliding DWs 
that was initially, for t — > —00, described with the field 
ingredient mi\ m x , is finally, for t — > 00, described 



with the field ingredient to 



+ ' 



Hence, during the 



7 



a) 



H 



-e- 



1©- 



4 



-®- 



©- 



-©- 



I 

-©- 



10- 



<8> 







8r 



-0* 



b) 



-0- 



-©- 



^0- 



-®- 



-0- 



-0- 



1©- 



.0' 







-0 1 



-0 1 - 



FIG. 4. Collision of DWs forming a hard magnetization bub- 
ble in a stripe. Scheme of the magnetization layout in the 
initial (final) state is shown in the left (right) picture for; a) 
pairs of transverse and vortex DWs of like chiralities and po- 
larities, b) pairs of transverse and vortex DWs of opposite 
chiralities and polarities. 



collision, DWs exchange their parameters xqi O X02 and 
their phase factors. The DWs reflect in a way that one 
can say they pass through each other without changing 
their velocities and polarities, however, with changing 
their character from the head-to-head one into the tail- 
to-tail one and vice versa and their chiralities (the process 
is illustrated in Fig. 4). This prediction corresponds to 
the result of the field-induced collision of a Bloch DW 
with a Neel DW in ID ferromagnet, [HI, as well as of the 
collision of spontaneously propagating (in absence of the 
dissipation) topological solitons (DWs) in ferromagnets 
which are intermediate structures between the Bloch and 
Neel DWs [13 . 

I have considered systems of infinite domains whose 
energy cannot be defined, however, the smaller a domain 
is the bigger percentage of the Zeeman part of its energy 
is lost per time unit due to the DW motion. The field- 
induced DW reflection induces the motion which contra- 
dicts the rule of energy minimization. Such a motion 
has to be decelerated and, eventually, suppressed when 
the decrease of the DW interaction energy equals the in- 
crease of the Zeeman energy. It results in formation of a 
stationary bubble which is a counterpart of hard bubbles 
in ID systems IS]. 



V. CONCLUSIONS 

In terms of the application of multi-DW systems to 
the magnetic storage, the important result of the present 
study is the prediction of bound states of noninteracting 
DWs [which are of opposite (like) chiralities and polar- 
ities] in the absence of any external field. Such bound 
states are stable with respect to simultaneous change of 
both the chirality and polarity of one of the DWs of the 



bubble. Unlike stationary bubbles in ID ferromagnets 
which are composed of one Neel DW and one Bloch DW, 
the bubbles in magnetic stripes are composed of DWs 
of the same energy, hence, they can be stable. In ID 
systems, the instability is a simple consequence of the 
fact that the Neel DWs are of higher energy than the 
Bloch ones and they tend to the reorientation into the 
Bloch walls in presence of fluctuations [36| . The possibil- 
ity of maintaining stable train of many DWs in magnetic 
stripes without external power supply (without applica- 
tion of the magnetic field) makes such systems potentially 
useful as magnetic information registers. 

The result of the field-induced collision of DWs of op- 
posite chiralities and polarities (or of like chiralities and 
polarities) is found to be their reflection. Hence, they can 
form a stripe counterpart of the hard bubbles of magne- 
tization of wide ferromagnetic platelets. 

Another finding to be stressed is the existence of sta- 
tionary bound state of two vortex DWs of like chiralities 
and of opposite polarities (a 2-7T-DW) in the absence of 
any external field. It follows from Fig. 3b that the en- 
ergy of the relevant pair of vortex DWs Eq(o) achieves 
minims at a =/= 0. They correspond to the final state of 
the long-term evolution of a breather. Such a state has 
no counterpart in ID ferromagnet while it has a counter- 
part in critical media described with Ginzburg-Landau 
(or a nonlinear Schrodinger) equation [l4| . 
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Appendix A: Estimation of boundary energy 

The magnetostatic energy of a magnetic element con- 
tains contributions that relate to interactions of surface 
charges, volume charges, and interaction between surface 
and volume ones 



E 



MS 



$^dS(x)dS(x 



cr(x)p(x 



rf dS(x)dV(x ), (Al) 



where p = — V • m, a = n • m. Following Ref. [371 ] . 
reducing one of the spatial dimensions with relevance to 
flat systems of thickness r and neglecting volume and 
base-surface terms, the above expression is transformed 
into the energy of the boundary of 2D system (up to the 
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multiplayer r) 



timate the boundary energy with 



tE b = t z I / cr(x)cr(x )ln(|x - x |/r)d/(x)dZ(x 
/ as J as 



cr(x)p(x ) + p(x)cr(x ) 



laS JSbas 

xln(|x-x'|/r)dZ(x)dS'(x'). 



(A2) 



Here Sbase denotes the surface of the platelet base. 
For any DW in the stripe, p(x) ~ —dm x /dx ~ 
[M 2 - m 2 x ] /(MS) with S denoting the DW width. I es- 



tEb=t z / / (n -m)(x)(n • m)(x )ln(|x — x |/r) 
'as J as 

>ii ' , , dm 



xd/(x)dZ(x ) 



'as Js base 
xln(|x-x'|/r)d/(x)d5(x') 

2<5r 2 ln((5/T) / m 2 (ir, 0, 0)dx 



(n • m)(x) 



dx 



x) 



-2at<;r 2 ln((5/i 



[Af 2 -m 2 (x,0,0)]d:c, (A3) 



where aw corresponds to an effective thickness of the 
surface layer of the stripe edge over which the magneti- 
zation is independent on normal coordinate z, (a <C 1 
and a oc 6). Since 5 ~ w, one arrives at Eb of ([5]). 



Appendix B: Explicit form of integrals 

/■°° / 

I ± (a, 9) = ( {sech 2 (y)tanh(— ?/ + a) — sech 2 (— y + a)tanh(y) =p cos(6*a) [— tanh(y) + tanh(— y + a)] 

J — DO 

xsech(y)sech(— y + a)} 2 + {— sech(y) [tanh(y)tanh(— y + a) + scch 2 (— y + a)] 

± cos(6'a)sech(— y + a) [tanh(j/)tanh(— y + a) + sech 2 (y)] } (Bl) 
+ {sin(6'a)sech(— y + a) [tanh(y)tanh(— y + a) + sech 2 (y)] } 2 

+ [sin(6>a)tanh(?/)sech(— y + a)] 2 + [sech(y)tanh(— y + a) ± cos(0a)tanh(y)sech(— y + a)] 2 ] dy 
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